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Abstract. 

We calculate the Casimir energy at spherical cavities within a host made up 
of an arbitrary material described by a possibly dispersive and lossy dielectric 
response. To that end, we add to the coherent optical response a contribution 
that takes account of the incoherent radiation emitted by the host in order to 
guarantee the detailed balance required to keep the system at thermodynamic 
equilibrium in the presence of dissipation. The resulting boundary conditions 
allow a conventional quantum mechanical treatment of the radiation within the 
cavity from which we obtain the contribution of the cavity walls to the density of 
states, and from it, the thermodynamic properties of the system. The contribution 
of the cavity to the energy diverges as it incorporates the interaction energy 
between neighbor atoms in a continuum description. The change in the energy 
of an atom situated at the center of the cavity due to its interaction with the 
fluctuating cavity field is however finite. We evaluate the latter for a simple case. 



1. Introduction 

Motivated by his finding tliat zero-point fluctuations may induce an attractive force 
between parallel conducting plates [T], Casimir proposed in 1956 that the zero-point 
force could be the Poincare stress involved a semiclassical model of the electron [2]. 
In that model the electron was considered as a spherical charge distribution stabilized 
by vacuum fluctuations. However, T. H. Boyer ^ showed in 1968 that the stress 
for a spherical conducting shell of radius a is indeed repulsive, since the Casimir 
energy turns out to be positive: E — .09235/2a. Subsequent calculations based on 
the Green's function method @], or a multiple scattering formalism confirmed 
Boyer's calculation. The more general problem of the Casimir effect of a dielectric 
ball was first considered by Milton [6] in absence of dispersion, and later on by 
Candelas [7] . Candelas argued that in presence of boundaries vacuum energies depend 
on a cutoff on transverse momenta, independently of the dielectric properties of the 
media. Therefore, Boyer's result should be corrected by terms associated to surface 
and curvature tensions. The Casimir forces for a dilute dielectric and diamagnetic 
sphere was studied by Brevik and Kolventvedt [8], and more recently by Klich [9], 
while the role of dispersion in this problem was discussed in Refs.|10| [TT| [12] . In 
this case, the Casimir stress may be attractive, but very sensitive to the specific 
values of the parameters characterizing the electric and magnetic response of the 
materials. An excellent review of different approaches to the Casimir effect of spherical 
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regions, together with appUcations to QCD bag models, higher dimensional spaces, or 
sonoluminescence can be found in Ref. 

In a series of papers jl41 1151 116] an expression for the Casimir force within planar 
cavities was derived without making particular models or assumptions about the 
nature of the walls. By considering that the system is in thermodynamic equilibrium, 
we obtained the energy and the stress tensor in a closed ancillary system that 
has the same optical response as the original system. This approach consistently 
incorporates evanescent fields and allows a fully quantum-mechanical treatment of the 
electromagnetic degrees of freedom. Unlike the calculations presented in [T31 [HI [TB] 
in Ref. 1 17] the fictitious system was eliminated, keeping only its essential property: 
that detailed balance should hold in thermodynamic equilibrium: for each photon 
that is not coherently reflected at the cavity walls and is therefore either absorbed or 
transmitted beyond the system, an identical photon has to be incoherently injected 
back into the cavity with no phase relation with the lost photon. In this paper we show 
that that approach may be straightforwardly generalized to study the Casimir effect 
of spherical cavities with arbitrary dielectric properties. We then apply the formalism 
to calculate the energy shift for a polarizable atom placed at the center of a cavity. 



2. Theory 

Consider a system made up of an empty spherical cavity of radius R carved out of 
an arbitrary material and with a scatterer situated at its center (Fig. [1]). Within the 
empty space of the cavity there are outgoing (o) and ingoing (i) electromagnetic waves 
with transverse electric (TE) and transverse magnetic (TM) polarization, described 
by the field 

Ff^ir-) = ff^hfikr)n,,{n), (1) 

where d = o,i describes the propagation direction, C — TE, TM describes the 
polarization, I — 0,1,2, . . . denotes the angular momentum, m = — L . . Z its projection 
along the z axis, h" = and hi = are the outgoing and the ingoing spherical 
Hankel functions respectively, k — lu/c is the wavenumber within vacuum, and we 

d TE "'d TE 

choose the scalar field = ■ f to describe the TE electromagnetic field 

Ef^, Bf^ and F;^™ = 4m™ ' ^ describe the TM electromagnetic field E™ , 
Bfjf [H]. When the outgoing radiation reaches the boundary of the cavity it is 
partially scattered back into the cavity with an amplitude and when the ingoing 

radiation hits the scatterer at the center it is scattered back towards the cavity with 
an amplitude that is, 

r = s^r, (2) 
r = scf, (3) 

where we removed the indices I, m, and C to simplify our notation. From Eqs. ([2]) 
and ([3]) we immediately obtain the normal modes of the system, given by 

1 - SfcSc = 0. (4) 

Taking into account the frequency {uj) dependence of Sf, and s^ we may solve Eq. Q 
to obtain the frequency spectrum w„ for each value of I, m, C. However, in the presence 
of absorbing materials or even of transparent, leaky materials, ujn would turn out to 
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Figure 1. Spherical cavity within arbitrary system with a scatterer at its center. 



be complex. The system would be open and ordinary quantum mechanics would not 
be applicable to the radiation field, i.e., TitUn would not be an energy quantum. 

The coefficients Sb and Sc describe the amplitude and the phase of the radiation 
that is coherently scattered back into the cavity. The energy that is not coherently 
scattered back into the cavity, described by 1 — |scP and 1 — |stp, is absorbed or 
it leaves the system through its external boundary. Nevertheless, in thermodynamic 
equilibrium, an absorbing scatterer or an absorbing enclosure has to eventually radiate 
back any radiation that it absorbs. The system should also admit photons from the 
vacuum that surrounds it to replenish those photons that were transmitted away. 
Detailed balance must hold and in the average, for each photon with numbers I, m, ( 
that leaves the cavity, an identical photon, with the same numbers but having no 
phase relation with the original photon, has to be injected back into the cavity. We 
mimic this incoherent radiation by a coherent field that is delayed a large time 
(x = b,c). To avoid interference with the coherently scattered field we take the limit 
— » oo. Thus, the incoherent radiation may be taken into account by replacing the 
scattering coefficients Sb and Sc by total scattering coefficients 

'-^^- = TTb^^ ^ = ^""^ 

Here, exp{iujT-)^) is the phase acquired during the large delay and is an extremely 
fast varying function of the frequency co, so that for any finite bandwidth interference 
effects would disappear. We assume that and are relatively slowly varying 
functions of frequency which are to be determined. The term in the numerator 
corresponds to the coherent scattering. The term a^^e"^-^" corresponds to re-radiation 
by the central scatterer (x = c), re-radiation by the walls of the cavity (x = b), or to 
photons that enter the system from outside to replenish the cavity losses. Finally, it 
may happen that a re-radiated photon fails to reach the cavity on its first attempt, as 
it may be re-absorbed or scattered away. Thus, we should allow for multiple injection 
attempts. These are accounted for by the term h^e^'^^^ in the denominator. 

As all the energy that leaves the cavity has to enter again in an equilibrium 
situation, the total scattering amplitudes must obey 

\S^? = 1, (6) 

which yields 



P + \a^\^ + 2'Res*^a^e^'^^ = 1 + l^^l^ + 2Reb^e''^^^. (7) 
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Figure 2. Integration contour 7 employed to obtain tlie density of states. We 
indicate the normal modes (crosses) separated approximately by Atii". 



Separating the slowly from the rapidly varying terms 

|sxl' + l«xl' = l + l&xl^ Re^flxe*"^- =Re6xe*"^-, (8) 



we obtain 

a^=e^^-, 6, = 4e'*-, (9) 

where are a slowly varying phases. The normal modes of the system in equilibrium 
are not given by Eq. (|4]) but by 

D = l- SbSc = 0, (10) 

which may be recast as 

2arg (1 + Sbe'^'e"^^") + 2arg (l + Sce'^^ e''^'^") 

+ {Sb + Sc)+u;{n + T,) = 2TTn, (11) 

with integer n. The first two terms in the LHS of Eq. (fTT|) oscillate rapidly, but are 
bounded within the interval (— 7r,7r), while the second term is slowly varying. Thus, 
the separation between nearby modes ojn is close to Aw" = 2tt/{Ti, + Tc) and the 
density of modes diverges in the limit Tb,Tc — *■ 00. This is to be expected, as our 
delayed re-radiation accounts implicitly for the interaction with a thermal bath which 
has infinite degrees of freedom. Our dissipative system together with the thermal 
bath forms an extended closed system whose modes are actually real fTUl and form a 
quasi-continuum. Our approach above is an alternative to the introduction of ancillary 
systems to account for dissipation [20l [16] . 

The actual number of modes AA'^ within a small frequency range Ail may be 
obtained using Cauchy's argument principle 

AN=^ i ^logfito), (12) 

where 7 is a contour that encircles counterclockwise the interval Af2, and /(cj) is an 
analytical function that has the same zeroes as D (Eq. pO)) ) and no poles within 
7. We choose a contour 7 that moves towards the right a distance ry below the real 
axis, then crosses the axis, moves back a distance rj above the real axis, and finally 
crosses the axis to closes upon itself |19J (Fig. [5]). Choosing as f{uj) the analytical 
continuation from the real axis unto the complex plane of 

f{uj) = (1 + S6e**^e*"^^)(l + s^e'^'^e''^'^'') 

- {sb + e'^-e"^^-)isc + e'^-e"^^"), (13) 
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we obtain 



AN 



1 
1 

+ ioj{Th 



A{log 



LU—ir) 

1 - 



+ i{&b + Sc) 



in the limit ^ oo, 77 — > 0, rfT^ — > c». Subtracting the number AA^o of modes 



corresponding to vacuum and the thermal reservoir only, obtained from Eq. p4[) by 
replacing ^ 0, Sc ^ 1, we obtain 

AiV - ATVo = --AImlog(l - SfesJ = p(w)A17, (14) 

TT 

where we identify the contribution of the scatterers to the density of states, 

p(u;) = -ilm-^log(l - SfeSc), (15) 



for each value of I, m, C- 

From the density of states, one can proceed to calculate all the thermodynamic 
quantities. For example, multiplying p{u)) by the ground state energy h(jj/2 of an 
oscillator of frequency co, integrating over all frequencies and adding over all angular 
momenta and polarizations we obtain the contribution of the scatterers to the ground 
state energy. 



l,mX 



du log 



1 - s 



(16) 



where we have already rotated the integration trajectory unto the imaginary axis. 

Eq. (|16[) and similar equations easily derived for other thermodynamic quantities 
are the main results of the present paper. In order to evaluate them the only 
requirement is knowledge of the scattering amplitudes corresponding to the surface of 
the cavity and to the scatterer at the center. 



3. Empty cavity within a uniform medium 

For an empty cavity = 1, as the incoming field becomes an outgoing field after 
crossing the origin. If the cavity is surrounded by a uniform medium with a dielectric 
function £(0;), then may be easily found by writing the field within the cavity as a 
linear combination of outgoing and ingoing fields (Eq. ([T])), writing the field within the 
medium as an outgoing field and matching both solutions through the usual boundary 
conditions, i.e., the continuity of the projections of both the electric and magnetic 
field along the surface. The result [H] is simply given by 

^b,lm\'^} QMoVi _ QViMo ' ^ 

where 

QAdA'd' ^ hf{kAR)Dhf{kA'R), A, A' = V,M, d,d' = i,o, (18) 
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kA is the wavenumber within vacuum {A = V , ky = k = co/c) or within the medium 
{A = M, kM = ^\/e), hf are the ingoing {d ~ i) or outgoing {d — o) spherical Hankel 
functions, and D is the operator 

bg{x) = g'{x)+glx. (19) 

Similarly, for TM polarization wc have 

^b,lm\^) — qMoVi _ QViMo ■ V^"'' 

As s^ is independent of m, we may replace the sum over m in Eq. (|16p by a factor 
of 21 + 1. 

In Fig. [3] we show the TM contribution to the integrand of Eq. as a function 
of the imaginary part of the frequency u = Lo/i calculated for a dielectric cavity 
with a Lorentzian response e = 1 + ujp/{ujQ — uP' — i^uj). The figure shows integrable 
singularities at w = 0. A second singularity is seen for even values of /. Nevertheless, 
it is clearly seen that the contributions to the energy grow with I and that Eq. (fT6|) 
doesn't converge. 

Our calculation above includes a realistic dielectric response for the medium 
surrounding the cavity and in our calculation the electromagnetic field permeates 
all space. Thus, the singularities above are of a different physical nature than 
the singularities arising in naive calculations for flat surfaces. For example, the 
singularities in the Casimir force among perfect conducting slabs may be removed by 
introducing a high frequency cutoff which accounts for the high frequency transparency 
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of real metals and when the mechanical properties of the field beyond the slabs is 
accounted for. Those recipes wouldn't cure the present divergence. 

Notice that large Vs correspond to spatial oscillations around the sphere with a 
small lengthscale d = 27ri?//. As we may assume a smallest lengthscale c?at of atomic 
dimensions, it seems reasonable to impose a corresponding cutoff at ^max — 27ri?/(iat. 
A careful analysis pi] shows that in this case Eq. (|16p converges, but that the leading 
terms are of order {R/datf and {R/datf', i.e., of the order of the number of atoms 
within the volume and the surface of a sphere of radius R. Mathematically, these 
terms arise from the leading terms in an expansion of the integrand of Eq. (|16p 
for large I, which, after summing over m are of order P and I (with logarithmic 
corrections). Physically, the reason is that Eq. (|16p includes the electromagnetic 
interaction between nearest neighbor atoms, which in a continuum description are 
infinitely small and infinitely close to each other. Similar divergent terms were found 
in a pairwise perturbative calculation for cavities surrounded by a very diluted 
dielectric. As argued in Ref. [S^, these short range contributions should be taken 
into account before conclusions about the sign of the Casimir force can be drawn. 

A finite Casimir energy may be obtained by subtracting from our result above 
those terms that contribute to the divergence before we take the limit dat 0. 
Nevertheless, the resulting energy would be an experimentally inaccessible quantity 
for a spherical cavity, as changing the cavity radius would require adding or 
removing atoms or else introducing strains that would produce an elastic stress [2 3 1 . 
Furthermore, a full calculation of the interatomic interactions would be required in 
order to compare experimental results to theoretical calculations, and our Eq. (fTH]) 
would be insufficient. One way out of these difficulties is to study other geometries 
where motion introduces no elastic stresses. For example, it has been found that 
for a piston sliding along a cylinder, the contribution of its position to the Casimir 
energy is free of cutoff dependent singularities and is attractive. 



4. Atom v^rithin a cavity 



On the other hand, Eq. pB]) would still be useful in situations where the geometry of 
the cavity is left unchanged. For example, it may be used to calculate the change 



27r ^ JO 



du log 



1 



c^Lm 



(iu) 



1 



(iu) 



(21) 



in the energy of the system when an atom is introduced at the center of the cavity. 
Using the usual selection rules, we obtain that the dispersion amplitude for an atom, 
^cim ~ 1' would be the same as for empty space, unless I = 1 and the polarization 
( = TM, in which case. 



TM 



lik^a 



1 



(22) 



where a is the electric-dipole polarizability of the atom. Therefore, there is only one 
finite term in Eq. (pij) . 

In Fig. [4] we show the contribution to the energy of a cavity from an atom lying 
at its center. The atomic polarizability is taken as a Lorentzian 

e^ /m 



a{Lj) 



(23) 
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Figure 4. Contribution to tlic energy of a polarizable atom lying at the center 
of a perfectly conducting cavity as a function of the radius R of the cavity. The 
polajrizability is taken as a Lorentzian with a resonance frequency ojo ■ 



with a single resonance frequency luq. The energy is negative and proportional to i? ^ 
and is of the order of tens of meV's for radii of a few nanometers. 

5. Conclusions 

Based on a scattering approach we have derived an expression for the Casimir energy 
of a spherical cavity with dispersive and absorptive dielectric properties. It turns 
out that the expression is divergent due to the summation over angular momenta, 

independently of the dielectric behavior of the cavity walls. This may be regularized 
by imposing a cutoff associated with the finite separation of atomic scatterers forming 
the boundary of the cavity, leading to contributions proportional to the number of 
atoms in the volume and surface of the cavity. On the other hand, the energy shift for 
an atomic radiator placed in the center of the cavity is finite, since it does not require 
to perform virtual work to modify the geometry of the cavity 
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